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^vq Abstract: We present an up-to-date global analysis of solar, atmospheric, reactor, and 

accelerator neutrino data in the framework of three-neutrino oscillations. We provide 
results on the determination of #13 from global data and discuss the dependence on the 

^ choice of reactor fluxes. We study in detail the statistical significance of a possible deviation 

of 023 from maximal mixing, the determination of its octant, the ordering of the mass states, 
and the sensitivity to the CP violating phase, and discuss the role of various complementary 
data sets in those respects. 
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1. Introduction 

It is now an established fact that neutrinos are massive and leptonic flavors are not sym- 
metries of Nature [1,2]. In the last decade this picture has become fully proved thanks to 
the upcoming of a set of precise experiments. In particular, the results obtained with solar 
and atmospheric neutrinos have been confirmed in experiments using terrestrial beams: 
neutrinos produced in nuclear reactors and accelerators facilities have been detected at 
distances of the order of hundreds of kilometers [3]. The minimum joint description of all 
these data requires mixing among all the three known neutrinos (z/ e , z/^, v T ), which can be 
expressed as quantum superpositions of three massive states fi (i = 1,2,3) with masses 
mj. This implies the presence of a leptonic mixing matrix in the weak charged current 
interactions [4, 5] which can be parametrized as: 



U 



' C12C13 S12C13 S 13 e~ l5cp 

-S12C23 - Cl2Sl 3 S23e*' 5cp Cl 2 C 2 3 - Sl2Sl 3 S23e i<5cp Cl 3 S 2 3 (1-1) 
\ Sl 2 S 2 3 - Ci2Si 3 C2 3 e j5cp -C12S23 - Sl2Sl 3 C 2 3e i<5cp C13C23 



where Cij = cos Oij and = sin0jj. In addition to the Dirac-type phase 5cp, analogous 
to that of the quark sector, there are two physical phases associated to the Majorana 
character of neutrinos, which however are not relevant for neutrino oscillations [6, 7] and 
are therefore omitted in the present work. Given the observed hierarchy between the solar 
and atmospheric mass-squared splittings there are two possible non-equivalent orderings 
for the mass eigenvalues, which are conventionally chosen as 

Am^i < (Ajn^ ~ Am^ > 0) ; (1.2) 
Amli < -(Amli ~ Am| 2 < 0) , (1.3) 
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with Am| = mf — m2 j- As it is customary we refer to the first option, Eq. (1.2), as Normal 
ordering (NO), and to the second one, Eq. (1.3), as Inverted ordering (10); in this form 
they correspond to the two possible choices of the sign of Am^. 1 In this convention the 
angles Oij can be taken without loss of generality to lie in the first quadrant, 0y G [0, vr/2], 
and the CP phase 5cp £ [0, 2-7r]. 

Within this context, Am^, [Am 2 ^!, 9\2, and #23 are relatively well determined, whereas 
till this year only an upper bound was derived for the mixing angle #13 and barely nothing 
was known on the CP phase <5cp and on the sign of Am^. This situation has dramatically 
changed with the data from the reactor experiments Daya Bay [8], Reno [9], and Dou- 
ble Chooz [10], which together with the increased statistics of long-baseline experiments 
T2K [11] and MINOS [12] have provided a a clear determination of the last unknown mixing 
angle #13. With these results at hand a first-order picture of the three-flavour lepton mixing 
matrix has emerged. Its precise determination, as well as that of the mass differences, can 
only be made by statistically combining the results of the oscillation searches. 

In this article, we present an up-to-date global analysis of solar, atmospheric, reactor 
and accelerator neutrino data in the framework of three-neutrino oscillations with the 
following outline. Sec. 2 contains the results of the global analysis and the extracted ranges 
of the oscillation parameters and the corresponding mixing matrix. In Sec. 3 we study 
the present determination of #13 and in particular we discuss the uncertainty associated 
to the choice of reactor fluxes. Sec. 4 focuses on our knowledge of 6*23 and <5cp with 
special emphasis on the present limitations of the statistical significance of the possible 
deviation of #23 from maximal, the determination of its octant and the sensitivity to the 
CP violating phase. In order to address these issues we study the role played by a potential 
complementarity of long baseline accelerator and reactor experiments (Sec. 4.1) and by the 
atmospheric neutrino results (Sec. 4.2). Finally, in Sec. 5 we summarize and discuss our 
results. Future updates of this analysis will be provided at the website Ref. [13]. Alternative 
recent global fits have been presented in Refs. [14, 15], see also [16, 17]. 

2. Oscillation Parameters: Results of the Global Analysis 

We include in our global analysis the results from Super-Kamiokande atmospheric neu- 
trino data from phases SK1-4 [18] (with addition of the 1097 days of phase SK4 over 
their published results on phases SK1-3 [19]). For what concerns long baseline acceler- 
ator experiments (LBL), we combine the energy distribution of disappearance events 
from K2K [20] with that obtained by MINOS in both [p v ) disappearance and v e (v e ) 
appearance with 10.8 (3.36) x 10 20 protons on target (pot) [21] (which update their pub- 
lished results [12,22]), and T2K v e appearance and disappearance data with phases 1-3, 
3.01 x 10 20 pot [23] (a factor ~ 2 increase with respect to their published results [11]), and 
phases 1-2, 1.43 x 10 20 pot [24,25], respectively. 

For oscillation signals at reactor experiments we include data from the finalized ex- 
periments CHOOZ [26] (energy spectrum data) and Palo Verde [27] (total rate) together 

1 In the following we adopt the (arbitrary) convention of reporting results for Am^ for NO and Am| 2 
for IO, i.e., we always use the one which has the larger absolute value. 
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Table 1: Three-flavour oscillation parameters from our fit to global data after the Neutrino 2012 
conference. For "Free Fluxes ± RSBL" reactor fluxes have been left free in the fit and short baseline 
reactor data (RSBL) with L < 100 m are included; for "Huber Fluxes, no RSBL" the flux prediction 
from [42] are adopted and RSBL data are not used in the fit. 



with the recent spectrum from Double Chooz with 227.9 days live time [28,29], and the 
total even rates in the near and far detectors in Daya Bay [30] with 126 live days of data (a 
factor 3 increase over their published results [8]) and Reno with 229 days of data-taking [9]. 
We also include the observed energy spectrum in KamLAND data sets DS-1 and DS-2 [31] 
with a total exposure of 3.49 x 10 32 target-proton-year (2135 days). 

Finally in the analysis of solar neutrino experiments we include the total rates from 
the radiochemical experiments Chlorine [32], Gallex/GNO [33] and SAGE [34]. For real- 
time experiments we include the 44 data points of the electron scattering (ES) Super- 
Kamiokande phase I (SKI) energy-zenith spectrum [35] and the data from the three phases 
of SNO [36-38] , including the results on the low energy threshold analysis of the combined 
SNO phases I— III [39]. We also include the main set of the 740.7 days of Borexino data [40] 
as well as their high-energy spectrum from 246 live days [41]. 

The results of the global analysis are shown in Figs. 1 and 2 where we show different 
projections of the allowed six-dimensional parameter space. The results are shown for two 
choices of the reactor fluxes as we will describe in more detail in the next section. The 
best fit values and the derived ranges for the six parameters at the la (3a) level are given 
in Tab. 1. For each parameter the ranges are obtained after marginalizing with respect to 
the other parameters. For sin 2 623 the la ranges are formed by two disconnected intervals 
in which the first one contains the absolute minimum and the second-one the secondary 
local minimum. Note that we marginalize also over the type of the neutrino mass ordering 
and the two local minima in sin 2 623 may correspond to different orderings. As visible in 
Fig. 2, for "Free Fluxes + RSBL" the best fit with sin 2 6*23 < 0.5 is for NO and the local 
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Figure 1: Global 3v oscillation analysis. Each panels shows two-dimensional projection of the 
allowed six-dimensional region after marginalization with respect to the undisplayed parameters. 
The different contours correspond to the two-dimensional allowed regions at la, 90%, 2a, 99% 
and 3tr CL (2 dof). Results for different assumptions concerning the analysis of data from reactor 
experiments are shown: full regions correspond to analysis with the normalization of reactor fluxes 
left free and data from short-baseline (less than 100 m) reactor experiments are included. For 
void regions short-baseline reactor data are not included but reactor fluxes as predicted in [42] are 
assumed. Note that as atmospheric mass-squared splitting we use Am^ for NO and Am| 2 for 10. 
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Figure 2: Global iv oscillation analysis. The red (blue) curves are for Normal (Inverted) Ordering. 
Results for different assumptions concerning the analysis of data from reactor experiments are 
shown: for solid curves the normalization of reactor fluxes is left free and data from short-baseline 
(less than 100 m) reactor experiments are included. For dashed curves short-baseline data are not 
included but reactor fluxes as predicted in [42] are assumed. Note that as atmospheric mass-squared 



splitting we use Am^ for NO and Am§ 2 for IO. 
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minimum with sin 2 #23 > 0.5 is for 10, while for "Huber fluxes, no RSBL" both minima 
are for 10. However, as also visible from the figure, the best fit values of the minima and 
the allowed range of sin 2 #23 are not very sensitivite to the mass ordering. The 3a range is 
connected. For Am^ (Am 2 2 ) the allowed ranges are formed by two disconnected intervals 
which correspond to the two possible mass orderings. From these results we conclude that: 

1. The present global analysis disfavours #13 = with a A% 2 ~ 100. This is mostly 
driven by the new reactor data from Daya Bay, Reno and Double Chooz (see Fig. 5 
in the next section). 

2. An uncertainty on $13 at the level of la remains due to a tension between predicted 
reactor neutrino fluxes and data from reactor experiments with baselines less than 
100 m, as we will discuss in detail in Sec. 3. 

3. Non-maximal #23 is favoured at the level of ~ 2a (~ 1.5er) for Normal (Inverted) 
ordering for either choice of the reactor fluxes. We elaborate more on this issue in 
Sec. 4. 

4. The statistical significance of the preference of the fit for the first octant of #23 
is < 1.5a (< 0.9<t) for Normal (Inverted) ordering for either choice of the reactor 
fluxes. 

5. When the normalization of reactor fluxes is left free and data from short-baseline (less 
than 100 m) reactor experiments are included, the absolute best-fit occurs for Normal 
ordering but the statistical significance of the preference Normal versus Inverted is 



6. The best fit occurs for Inverted ordering when reactor short-baseline data are not 
included and reactor fluxes as predicted in [42] are assumed but the statistical sig- 
nificance of the preference Inverted versus Normal is < 0.75o\ 

7. The statistical significance of the effects associated with 5cp is < 1.5a (< 1.75cr) for 
Normal (Inverted) ordering (see Sec. 4). 

From the global x 2 an d following the procedure outlined in Ref. [43] one can derive 
the following 3a CL ranges on the magnitude of the elements of the leptonic mixing matrix 



By construction the derived limits in Eq. (2.1) are obtained under the assumption of the 
matrix U being unitary. In other words, the ranges in the different entries of the matrix 
are correlated due to the fact that, in general, the result of a given experiment restricts 
a combination of several entries of the matrix, as well as to the constraints imposed by 
unitarity. As a consequence choosing a specific value for one element further restricts the 
range of the others. 



< 0.7cr. 




0.513 
0.416 
0.409 



0.585 
0.730 
0.725 




(2.1) 
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3. Determination of 613 and Flux Uncertainties 

At present, reactor experiments with L ~ 1 km provide us with the most precise determina- 
tion of #13. Up to very recently the interpretation of neutrino oscillation searches at nuclear 
power plants was based on the calculations of the reactor v e flux from Ref. [44]. Indeed, 
the observed rates at all reactor experiments performed so-far at distances L < 1 km are 
consistent with these fluxes, therefore setting limits on v e disappearance. Over the last two 
years the flux of v e emitted from nuclear power plants has been re-evaluated [42,45], yield- 
ing roughly 3% higher neutrino fluxes than assumed previously. As discussed in Ref. [46] 
this might indicate an anomaly in reactor experiments at L < 1 km, which according to 
the new fluxes observe a slight deficit. For the Chooz and Palo Verde experiments at 
L ~ 1 km a non-zero #13 could lead to v e disappearance accounting for the reduction of 
the rate. However, Am 2 3 and #13 driven oscillations will have no effect in short-baseline 
(SBL) experiments with L < 100 m. 

Motivated by this situation we follow here the approach of Ref. [47] and study the 
dependence of the determined value of #13 on the assumptions about the reactor fluxes, see 
also [48]. In particular we consider the impact of allowing for a free normalization of the 
reactor fluxes / nux and/or of including in the analysis the results of the reactor experiments 
with L < 100 m Bugey4 [49], ROVN04 [50], Bugey3 [51], Krasnoyarsk [52,53], ILL [54], 
Gosgen [55], SRP [56], and ROVNO88 [57], to which we refer as reactor short-baseline 
experiments (RSBL). 

The outcome is illustrated in Fig. 3 where in the upper panels we show Ax 2 from the 
Chooz and Palo Verde, Double Chooz, Daya Bay and Reno experiments as a function of 
sin 2 #13 for various assumptions on the fluxes. As seen in the left upper panel, when the 
new fluxes are taken at face value and RSBL reactor experiments are not included in the fit, 
even Chooz and Palo Verde would prefer #13 > but as soon as RSBL reactor experiments 
are included in the fit, the preference essentially disappears [47] independently of whether 
the normalization of the fluxes is also left as a free parameter. The lower left panel shows 
the contours in the plane of sin 2 #13 and the flux normalization /fj ux for the analysis with 
free reactor flux normalization with and without including RSBL data. Central panels 
show the dependence of the value of #13 obtained from the analysis of Double Chooz. As 
seen in the figure, the best fit value as well as the statistical significance of the non-zero #13 
at Double Chooz severely depends on the reactor flux assumption. This is due to the lack 
of the near detector in Double Chooz at present. 2 Conversely the right panels show the 
results of Daya Bay and Reno which rely on a near-far detector comparison and therefore 
are independent of the overall flux normalization. 3 Once they are included in the global 
analysis the impact of RSBL data is reduced but it still makes a difference in the final best 
fit values and ranges. To quantify this effect we chose to show the results of the global 

2 Let us mention that the official Double Chooz analysis [10, 28] adopts the Bugey4 measurment of the 
reactor flux as normalization for the reactor flux, and therefore it is similar to our "Free Fluxes + RSBL" 
analysis. 

3 Note that to-date neither Daya Bay nor Reno have published results on the absolute flux determination. 
Therefore we always have to include free normalization factors in the fit to their data. 
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Figure 3: Upper: A% 2 as a function of sin #13 for the different reactor experiments and different 
assumptions on the fluxes as labeled in the figure. In this figure we fix Am^ = 2.47 x 10 -3 eV . 
Lower: contours in the plane of sin 2 # 13 and the flux normalization /a ux . Full regions (lines) 
correspond to analysis with (without) including the RSBL experiments. 



analysis for the two limiting assumptions of either taking the predicted fluxes (with the 
related uncertainties and correlations) of [42] and ignore the RSBL data (which we label in 
the figures as "Huber" ) or to allow for a free normalization of the reactor fluxes and include 
the RSBL data to reduce its possible allowed range (labeled as "Free Fluxes + RSBL" ) . 

It is also interesting to notice that since the dominant oscillation probability in these 
reactor experiments with L ~ 1 km is 

P VM = 1 - sin 2 29 13 sin 2 (^f^) + 0(a 2 ) , (3.1) 

with a = Am^/Amjj, then the rates observed in the detectors at different baselines can 

4 In the numerical analysis higher order effects associated to Am|i are included, and they have a notice- 
able effect on the extraction of #13, especially for Daya Bay. In principle those effects could even distinguish 
Normal and Inverted orderings [58-60], an effect which however is below the present sensitivity of the 
experiments. 
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Figure 4: 3er allowed regions in the plane of lAmjj-J and sin 2 # 13 for different combinations of the 
reactor experiments. The region labeled "ALL REACTORS" does not include Kamland. 

provide an independent determination of Am^ [61]. We show in Fig. 4 the 3<t allowed 
regions in the plane | Am 2 ^ versus sin 2 #13 for different combinations of the reactor exper- 
iments. Of course the accuracy on lAm^j from these data is much worse than the one 
from MINOS (although consistent). One may expect improved sensitivity for IA777.3-J from 
reactors once spectral data from Daya Bay and Reno will be included. 

We conclude this section by presenting in Fig. 5 the dependence of Ax 2 on sin 2 #13 
for the different data samples included in the global analysis. In the upper left panel 
we summarize the present status of the "hint" for a non-zero #13 from the mismatch of 
the best-fit point values of Am^ and #12 between the solar analysis and KamLAND in 
the framework of 2-v oscillations [62-67]. As discussed in detail in Ref. [67], and seen in 
the figure, the statistical significance of this effect depends on the Standard Solar Model 
employed in the analysis, either the GS98 model or the AGSS09 [68] 5 and, to lesser extend, 
on the capture cross-section in gallium employed 6 . As seen in the lower panel of Fig. 5, 

5 GS98 is based on the older solar abundances leading to high metallicity and which perfectly agreed 
with helioseismological data. AGSS09 uses the new precise determination of the solar abundances which 
imply a lower metallicity and cannot reproduce the helioseismological data. This conflict constitutes the 
so-called "solar composition problem". 

6 To explain the lower-than-expected rate observed in the calibration of the GALLEX and SAGE detectors 
(the so-called "Gallium anomaly" ) Ref. [34] considers a modification of the capture cross-section of Ref. [69] 
where the contribution from the two lowest-lying excited states in 71 Ge is set to zero. In our analysis 
labeled "GS98" we adopt the cross section from [69], whereas for the analysis labeled "AGSS09" we used 
the modified cross section following [34], see [67] for details. 
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Figure 5: Dependence of A\ 2 on sin 2 #13 for the different data samples and assumptions as labeled 
in the figure and the corresponding la ranges. 

once the solar and KamLAND results are combined with the oscillation signatures from 
reactor and LBL experiments the final allowed range of the oscillation parameters becomes 
robust under changes on these factors on the solar analysis. 



4. Determination of #23 and Sep 

As seen in Sec. 2 our results show that non-maximal 823 is favoured at the level of ~ 2a 
(~ 1.5a) for Normal (Inverted) ordering while the statistical significance of the effects 
associated with 5cp is < 1.5a (< 1.75c) for Normal (Inverted) ordering. Next we discuss 
the relative role of LBL, reactor and atmospheric neutrino data in these results. 

4.1 The beam reactor interplay 

Since the advent of data on — > v e searches from T2K and MINOS on the one side, 
and #13 reactor experiments Double Chooz, Daya Bay, and Reno on the other side, the 
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long anticipated complementarity of beam and reactor experiments [70,71] is now a reality. 
As shown in Sec. 2, the global analysis indicates a deviation of 023 from the maximal 
mixing value of 45°, roughly at the level of 1.7cr — 2a, see Fig. 2. If confirmed, such a 
deviation would have profound implications for neutrino mass models based on flavour 
symmetries. An important contribution to this effect comes from recent MINOS data on 
disappearance. Neglecting effects of Am^ and the matter effect, the relevant survival 
probability in MINOS is given by 

P v ^ = l-4|C/ At3 | 2 (l-|[/ M 3| 2 )sin 2 ^^, |^ 3 | 2 =sin 2 2 3cos 2 13 , (4.1) 

where L is the baseline and E is the neutrino energy. Hence, the probability is symmetric 
under It/^) 2 — > (1 — |£^3| 2 )- In the two-flavour limit of #13 = this implies that the data 
is sensitive only to sin 2 2023, which for #23 7^ 45° leads to a degeneracy between the first 
and second octants of #23 [72]. Indeed, recent data from MINOS [21] give a best fit point 
of sin 2 20 ~ 0.94 if analyzed in a two-flavour framework. 

Since #13 is large one can try to explore a synergy between long-baseline appearance 
experiments and an independent determination of #13 at reactor experiments in order to 
resolve the degeneracy [70,72,73]. Let us look at the appearance probability relevant for the 
v u. v e searches at T2K and MINOS. Expanding to second order in the small parameters 
sin 0i3 and a = Am^/Am^ and assuming a constant matter density one finds [74-76]: 



a ■ 2n ■ 2n sin 2 A(l - ,4) 2 • 2 o/i 2n sin 2 AA 
p v^u e ~4 sin 13 sin 23 — ^ _ ^ 2 V a sin 20 i2 cos 23 — ^ — 

r s sin AA sinA(l - A) 
+ 2 a sin i3 sin 20 i2 sin 20 23 cos(A ± 5 C p) 



(4.2) 



A I- A 



with the definitions 

A?n 2 L 2EV 



A==^JK A=j^-, (4.3) 
4ii Anijj 

where L is the baseline, E is the neutrino energy, and V is the effective matter potential [77]. 
Note that a, A, and A are sensitive to the sign of Am^ (i.e., the type of the neutrino 
mass ordering). The plus (minus) sign in Eq. (4.2) applies for neutrinos (antineutrinos), 
and for antineutrinos V — > —V, which implies A — > —A. It is clear from Eq. (4.2) that in 
the case of large matter effect, A > 1, the terms (1 — A) depend strongly on the type of 
the mass ordering, and for A = 1 (possible for neutrinos and NO, or anti-neutrinos and 
10) a resonance is encountered [78]. Numerically one finds for a typical matter density of 
3 g/cm 3 



14* 0-09 (Ai)( lAm \ l 2 ) ■ (4.4) 
1 1 V GeV / \2.5 x IO- 3 eV 2 / V ; 

Since for T2K E ~ 0.7 GeV, matter effects are of order few percent, whereas in experiments 
like NOvA [79] with E ~ 2 GeV we can have \A\ ~ 0.2. Note that a 2 fx 10~ 3 , which implies 
that the second term in the first line of Eq. (4.2) gives a very small contribution compared 
to the other terms. An important observation is that the first term in Eq. (4.2) (which 
dominates for large 013) depends on sin 2 023 and therefore is sensitive to the octant. Reactor 
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Figure 6: Left: Preferred regions at 68% CL in the sin 2 2# 13 — Sep plane. The contour curves 
correspond to T2K + MINOS appearance data, where sin 2 6*23 is fixed to the two degenerate solution 
in the 1st (red) and 2nd (blue) octant. We define contours for 2 dof with respect to the global 
minimum which is indicated by a triangle. The gray region corresponds to the #13 determination 
from the reactors Double Chooz, Daya Bay, Reno (la band for sin 2 13 , 1 dof). Right: x 2 (<5) from 
beams (dashed) and beams+reactors (solid) with the same color coding as in the left panels. The 
solid curves are computed by adding Axg la = (sin 2 6" 13 - 0.023) 2 /(0.0023) 2 to the x 2 horn T2K and 
MINOS appearance data. Upper (lower) panels are for NO (10). The other oscillation parameters 
are fixed to the best fit values from Tab. 1 (Free Fluxes + RSBL). 



experiments with L ~ 1 km, on the other hand, provide a measurement of #13 independent 
of 623, see Eq. (3.1). Hence, by combining the data from reactor experiments such as 
Double Chooz, Daya Bay, and Reno with the appearance data from T2K and MINOS one 
should be sensitive in principle to the octant of #23. 

The situation from present data is illustrated in Fig. 6, where we show the determina- 
tion of #13 from the beam experiments T2K and MINOS as a function of the CP phase <5cp 
and the octant of #23 , where we have chosen values motivated by the MINOS disappear- 
ance result. The resulting regions in sin 2 2#i3 are compared to the reactor measurements 
from Double Chooz, Daya Bay, and Reno. It is clear from that figure that for present data 
from beams and reactors it is not possible to distinguish between 1st and 2nd 623 octants. 
For both possibilities overlap regions between beams and reactors can be found, although 
at different values of 5qp- Therefore, current data from reactor and long-baseline beam 
experiments are not able to resolve the 623 octant degeneracy. The lifting of the degeneracy 
(at low CL) visible in Fig. 2 appears due to atmospheric neutrino data, to be discussed 
below. 
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In principle the reactor-beam combination should also offer some sensitivity to the 
CP phase 5qp- This is shown in the right panels of Fig. 6. We see that if the octant of 
023 and the neutrino mass ordering were known, already present data from the beam and 
reactor experiments used in that figure would show quite sizeable dependence on the CP 
phase, depending on which of the 4 degenerate solutions is considered. However, it is also 
clear from the figure that once we marginalize over those four solutions, x 2 (^cp) becomes 
very flat and essentially all values of <5cp would be consistent within A% 2 < 0.5. This is 
a real-life example of how degeneracies can seriously spoil the sensitivity of long-baseline 
data [80]. The somewhat larger 5cp dependence visible in Fig. 2 follows again from the 
global fit including atmospheric neutrinos, as discussed next. 

4.2 The impact of atmospheric neutrinos 

Atmospheric neutrinos provide a powerful tool to study neutrino oscillations. The neutrino 
source contains v e and Vn as well as neutrinos and antineutrinos and furthermore, for a 
good fraction of the events, neutrinos travel through the Earth matter. In the context of 
3is mixing, the dominant oscillation channel of atmospheric neutrinos is — > v T driven 
by [Am^-J with an amplitude controlled by 023- However, the richness of the atmospheric 
neutrino beams and baselines opens up the possibility of sensitivity to subleading oscillation 
modes, triggered by Am^ and/or #13, especially in the light of the new determination of 
#13. In particular, they can shed light on the octant of 6*23 and perhaps on the value of <5cp 
and the type of neutrino mass ordering. 

An interesting observable is the excess of e-like events (relative to the no-oscillation 
prediction N®), since in the two-flavour limit one expects N e = N® and therefore any 
deviation of the observed number of events from N® should be due to subleading effects. 
Such excess can be written in the following way (see, e.g., [81]): 

^ - 1 ~ ("in 2 23 - l)P 2 v(Aml 1 ,e 13 ) 

+ (rcos 2 e 23 -l)P 2l/ (Am 2 21 ,e 12 ) (4-5) 

- sin 6>i3 sin 20 23 r Re(A* ee A^ e ) . 

Here r = $^/$ e is the flux ratio with r ~ 2 in the sub-GeV range and r ~ 2.6 — > 4.5 
in the multi-GeV range. P 2v (A.m 2 ,9) is an effective two-flavour oscillation probability 
and A ee ,A^ e are elements of a transition amplitude matrix. The three terms appearing 
in Eq. (4.5) have a well defined physical interpretation. The first term is important in 
the multi-GeV range and is controlled by the mixing angle $13 in P 2u (Am 2 1 , 013). This 
probability can be strongly affected by resonant matter effects [82-86]. Depending on 
the mass hierarchy the resonance will occur either for neutrinos or antineutrinos. The 
second term is important for sub-GeV events and it takes into account the effect of "solar 
oscillations" due to Am|i and 0i2 [87-90]. Via the pre-factor containing the flux ratio 
r both, the first and second terms in Eq. (4.5) depend on the octant of 023, though in 
opposite directions: the multi-GeV (sub-GeV) excess is suppressed (enhanced) for 023 < 
45°. Finally, the last term in Eq. (4.5) is an interference term between 0i3 and Am^ 
amplitudes and this term shows also dependence on the CP phase 5qp [81,90]. 
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Subdominant three neutrino effects can also affect //-like events. For example for 
multi-GeV muon events one can write the excess in //-like events as [91,92] 



The first term is controlled by #13 and is subject to resonant matter effects, similar to the 
first term in Eq. (4.5), though with a different dependence on #23 and the flux ratio. In the 
second term, ^33 is a probability amplitude satisfying /^(Am^, #13) = 1 — |^33| 2 . In the 
limit #13 = we have Re(^33) = cos(Am| 1 L/2£'), such that the second term in Eq. (4.6) 
just describes two-flavour z/ M — > vacuum oscillations. 

The statistical significance of these effects in the present global analysis is shown in 
Fig. 7 where we show the dependence of A% 2 on sin 2 623 and 5qp- The three curves in each 
panel correspond to the global analysis including atmospheric data from phases SK1-4 
(red full lines), the analysis without including the atmospheric neutrino data (black full 
lines), and the global analysis with the previous atmospheric sample from phases SK1- 
3 (green dashed lines). For the sake of comparison with the results of the analysis of 
the Super-Kamiokande collaboration, in our atmospheric analysis with phases SK1-4 the 
no-oscillation expectations have been normalized to those obtained with the new Honda 
fluxes [93], while the analysis from SK1-3 is done with the previous set of fluxes from 
the same group [94]. For simplicity we only show the results of the analysis with the 
normalization of reactor fluxes is free and data from short-baseline reactor experiments 
included (Free Fluxes + RSBL). Similar behaviour is found for the analysis using reactor 
fluxes as predicted in [42] (Huber). 

The figure illustrates that in the analysis without atmospheric neutrinos (full black 
lines) a preference for non-maximal value of 623 at the level of 1.7a (2a) for Normal (In- 
verted) ordering is observed. Such result is mainly driven by the new MINOS disappear- 
ance data [21]. On the other hand, we do not observe any statistically relevant sensitivity 
to the octant of #23 from the analysis without atmospheric neutrinos. This is, as mentioned 
in the previous section, with the present data from beams and reactors it is not possible 
to distinguish between first and second #23 octants. 

Comparing black with either red or green lines in Fig. 7 we see that in the global 
analysis including atmospheric data a preference for the first octant is observed. This can 
be attributed to a zenith-angle independent event excess in the sub-GeV e-like data in 
Super-Kamiokande. Such an excess can be explained by oscillations due to Am^ [87-89]. 
For sub-GeV events the second term in Eq. (4.5) is relevant. In that energy regime r ~ 2 
and for sin 2 #23 ~ 0.5 the pre-factor (r cos 2 #23 — 1) is suppressed, whereas in the first octant 
with sin 2 #23 < 0.5 an excess is induced. We also see that despite the larger statistics, the 
effect is slightly less significant in the analysis including the SK4 sample (full red lines) 
versus the one without that sample (dashed green lines). This is in agreement with the 
preliminary analysis performed by the Super-Kamiokande [18] collaboration of their full 
phases SK1-4 in which this excess is less significant compared to phases SK1-3. Let us 




(4.6) 



- sin 2 2023 [1 - Re(4$3)] 
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Figure 7: Ax 2 as a function sin 2 #23 and Sqp for three different analyses as labeled in the figure. 
Upper (lower) panels correspond to Normal (Inverted) ordering. 

also mention that due to the way the data are presented we have to use different binning in 
fitting SK1-3 compared to SK1-4, in particular of the sub-GeV samples, which also affects 
the sensitivity to the 623 octant. 

Given the preference for the first #23 octant once atmospheric data is included, also 
the sensitivity to <5cp is somewhat increased as seen in the lower panels of Fig. 7. This can 
be understood from both the effect of 5qp in the atmospheric data, as well as in the long- 
baseline experiments. Looking at Fig. 6 we see that by once the solution with sin 2 623 < 45° 
is favoured, the beam-reactor combination provides a better sensitivity to Sqp, visible in 
the right panels. This results in the final sensitivity shown in Fig. 7, which is at the level of 
Ax 2 ~ 3. We emphasize again the crucial interplay of different data sets necessary for this 
sensitivity to emerge: MINOS v„ disappearance prefers sin 2 2#23 < 1, atmospheric data 
disfavours sin 2 Q23 > 0.5, and the atmospheric data together with the i>n — > v e data from 
beams combined with the #13 determination from reactors provide sensitivity to <5cp. In 
brief, in the present global analysis both the sensitivity to the octant of #23 and to <5cp still 
relies on the observability of sub-dominant oscillation effects in the atmospheric neutrino 
analysis. 

In this respect it is important to stress that already since SK2 the Super-Kamiokande 
collaboration has been presenting its experimental results in terms of a large number of 
data samples. The rates for some of those samples cannot be theoretically predicted (and 
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therefore included in a statistical analysis) without a detailed simulation of the detector, 
which can only be made by the experimental collaboration itself. Hence, although our 
results represent the most up-to-date analysis of the atmospheric neutrino data which can 
be performed outside the collaboration, such an analysis has unavoidable limitations. For 
example, as can be seen from Eqs. (4.5) and (4.6) there can be some features in e-like or 
/i-like data samples which exhibit a different dependence on #23, and which of those subtle 
effects dominates depends on details of the detector simulation, binning, and treatment of 
systematic uncertainties. 

5. Discussion 

We have presented the results of an updated analysis of solar, atmospheric, reactor and 
accelerator neutrino data in the framework of three-neutrino oscillations. Quantitatively 
the results of the present determination of the oscillation parameters is listed in Table 1. 
The corresponding leptonic mixing matrix is given in Eq. (2.1). 

At present the most important players in the new determination of #13 are the reactor 
experiments (see Fig. 5). The global fit excludes #13 being zero with Ax 2 ~ 100. The results 
of reactor experiments without a near detector, in particular Double Chooz, Chooz and Palo 
Verde, depend on the expected rates as computed with some prediction for the neutrino 
fluxes from the reactors. This brings up the possible uncertainty of the determination of #13 
associated with the tension between the recent reevaluation of the reactor neutrino fluxes 
in Refs. [42,45] and data from reactor experiments with baselines less than 100 m. In Sec. 3 
we have quantified this effect and conclude that an uncertainty on the determination of #13 
at the level of la remains due to this so-called reactor anomaly [46]. As the statistics of 
Daya Bay and Reno experiments increases and with the operation of the near detector in 
Double Chooz this contribution to the error is expected to decrease. 

Sec. 4 focuses on our present knowledge of #23 and 5qp with special emphasis on 
the present limitations of the statistical significance of the possible deviation of 623 from 
maximal, the determination of its octant and the sensitivity to the CP violating phase. 
In particular we have studied the role played in the present analysis by LBL, reactor 
and atmospheric neutrinos on these issues. We find that the global analysis prefers a non- 
maximal value of 023 at the level of 1.7-2<r and that this result is mostly driven by the recent 
MINOS disappearance data. However we do not observe any sensitivity to the octant 
of #23 from the global analysis without atmospheric neutrinos. Although the combination 
of LBL and reactor results has the potential to disentangle the octant of Q23 , within their 
present precision this effect is not statistically significant as illustrated in Fig. 6. The lifting 
of the octant degeneracy (which we find to be at most a 1.5cr effect) visible in Fig. 2 appears 
due to atmospheric neutrino data (see also Fig. 7). 

Equivalently, in principle the reactor-beam combination could also offer some sensi- 
tivity to the CP phase. However, we find that this effect is not significant enough and 
without inclusion of atmospheric neutrinos in the analysis all values of 5qp are consistent 

7 For details on our simulation of the data samples and the statistical analysis see the Appendix of 
Ref. [3]. 
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within Ax 2 < 0.5 (Fig. 7). The observed sensitivity to <5cp hi the global analysis is a com- 
bined effect of (i) MINOS disappearance favouring sin 2 2^23 < 1, (h) atmospheric data 
disfavouring sin 2 #23 > 0.5, and (hi) the — > v e data from LBL experiments combined 
with the #13 from reactors. Altogheter they provide a statistical significance for <5cp at the 
level of 1.7o". 

We conclude that within the present accuracy of the LBL experiments and reactor 
experiments, the small statistical significance of the determination of the octant of #23 
and of <5cp is driven by the observation of their subdominant effects on the atmospheric 
neutrino events. Finally we mention that, both neutrino mass orderings (Normal and 
Inverted) provide a fit of very similar quality to the global data, with A% 2 ~ 0.5. 

Future updates of this analysis will be provided at the website Ref. [13]. 
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